library(plotly)
data <- read.csv("TB_burden_countries_2019-12-23.csv")
data2018<-data[which(data$year==2018), ]
head(data2018)
# light grey boundaries
l <- list(color = toRGB("grey"), width = 0.5)
# specify map projection/options
g <- list(
showframe = FALSE,
showcoastlines = FALSE,
projection = list(type = 'Mercator')
)
p <- plot_geo(data2018) %>%
add_trace(
z = ~e_pop_num, color = ~e_pop_num, colorscale = 'Viridis',
text = ~country, locations = ~iso3, marker = list(line = l)
) %>%
colorbar(title = 'Population', tickprefix = '') %>%
layout(
title = '2019 Tuberculosis (TB)<br>Source:<a href="https://www.who.int/tb/country/data/download/en/">World Health Organization</a>',
geo = g
)
p
# light grey boundaries
l <- list(color = toRGB("grey"), width = 0.5)
# specify map projection/options
g <- list(
showframe = FALSE,
showcoastlines = FALSE,
projection = list(type = 'Mercator')
)
p <- plot_geo(data2018) %>%
add_trace(
z = ~e_mort_exc_tbhiv_100k_lo, color = ~e_mort_exc_tbhiv_100k_lo, colorscale = 'RdBu',
text = ~country, locations = ~iso3, marker = list(line = l)
) %>%
colorbar(title = 'Case Count', tickprefix = '') %>%
layout(
title = '2019 Tuberculosis (TB)<br>Source:<a href="https://www.who.int/tb/country/data/download/en/">World Health Organization</a>',
geo = g
)
p
# light grey boundaries
l <- list(color = toRGB("grey"), width = 0.5)
# specify map projection/options
g <- list(
showframe = FALSE,
showcoastlines = FALSE,
projection = list(type = 'Mercator')
)
p <- plot_geo(data2018) %>%
add_trace(
z = ~e_mort_exc_tbhiv_num, color = ~e_mort_exc_tbhiv_num, colorscale = 'Magma_r',
text = ~country, locations = ~iso3, marker = list(line = l)
) %>%
colorbar(title = 'Case Count', tickprefix = '') %>%
layout(
title = '2019 Tuberculosis (TB)<br>Source:<a href="https://www.who.int/tb/country/data/download/en/">World Health Organization</a>',
geo = g
)
p
# light grey boundaries
l <- list(color = toRGB("grey"), width = 0.5)
# specify map projection/options
g <- list(
showframe = FALSE,
showcoastlines = FALSE,
projection = list(type = 'Mercator')
)
p <- plot_geo(data2018) %>%
add_trace(
z = ~cfr_hi, color = ~cfr_hi, colorscale = 'YlOrRd',
text = ~country, locations = ~iso3, marker = list(line = l)
) %>%
colorbar(title = 'Casea Count', tickprefix = '') %>%
layout(
title = '2019 Tuberculosis (TB)<br>Source:<a href="https://www.who.int/tb/country/data/download/en/">World Health Organization</a>',
geo = g
)
p
# light grey boundaries
l <- list(color = toRGB("grey"), width = 0.5)
# specify map projection/options
g <- list(
showframe = FALSE,
showcoastlines = FALSE,
projection = list(type = 'Mercator')
)
p <- plot_geo(data2018) %>%
add_trace(
z = ~cfr_lo, color = ~cfr_lo, colorscale = 'Greens',
text = ~country, locations = ~iso3, marker = list(line = l)
) %>%
colorbar(title = 'Casea Count', tickprefix = '') %>%
layout(
title = '2019 Tuberculosis (TB)<br>Source:<a href="https://www.who.int/tb/country/data/download/en/">World Health Organization</a>',
geo = g
)
p